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Abstract 

In this work we calculate the ^5*0 gap energies of A hyperons in neutron star 
matter. The calculation is based on a solution of the BCS gap equation for 
an effective G-matrix parameterization of the A— A interaction with a nuclear 
matter background, presented recently by Lanskoy and Yamamoto. We find 
that a gap energy of a few tenths of MeV is expected for A Fermi momenta 
up to about 1.3 fm^^. Implications for neutron star matter are examined, 
and suggest the existence of a A ^5*0 superfluid between the threshold baryon 
density for A formation and the baryon density where the A fraction reaches 
15-20%. 



PACS: 26.60. +c, 97.60.Jd, 14.20. Jn 
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I. INTRODUCTION 



The theory of neutron star structure directly relates global properties of these stars to 
various aspects of many-baryon physics. One fundamental issue is whether pairing forces 
among the baryons can give rise to baryon superfluids in the inner crust and quantum 
cores of neutron stars. While nucleon pairing in neutron stars has received much attention, 
quantitative estimates of pairing of other baryon species has not been performed to date, due 
to lack of relevant experimental data. In this work we use some recent analysis of hypernuclei 
to make a first attempt at determining superfluid gaps for A hyperons in neutron star matter. 

Since first suggested by Migdal nucleon pairing in nuclear matter has been the subject 
of many studies. Both former and recent [4-8] works typically find ^5*0 neutron pairing 
for neutron matter density, p„, in the range of O.lpo < Pn < 0.5po, where po ~ 0.16 fm~^ 
is the nuclear saturation density. At higher densities, the ^5*0 interaction turns repulsive, 
and pairing is possible through higher order interactions, mainly [p|,p!0| . The energy gap 
found for the ^5*0 neutron superfluid is typically of the order of a few MeV, although recent 
works 1^,0] suggest that quasi-particle correlations could lower the energy to about 1 MeV. 
Estimates of the gap are typically of the order of a few tenths of a MeV. It should be noted 
that published results for the pairing energy gaps differ by as much as a factor of three. The 
difficulty in obtaining accurate results is mainly due the the problem of consistently including 
background medium effects. Uncertainties in the two-body interactions pose an additional 
problem. 

As the temperature of a neutron star is expected to drop below 0.1 MeV (~ 10^ K) 
within about one day from its birth, it is widely accepted that nucleon superfluids exist in 
different regions of the star. The qualitative picture of a neutron star includes a ^5*0 neutron 
superfluid in its inner crust (along with neutron rich nuclei), and a neutron superfluid 
in the quantum liquid core. The protons in the core, having a partial density of about 10% 
of the neutrons, are also expected to be in a ^S'o superfluid, with an energy gap of about 1 
MeV. m- 

Baryon superfluids are expected to have a number of important consequences on neutron 
star physics including several observational effects, such as pulsar glitch phenomena and 
cooling rates. The crustal neutron superfluid is expected to play a incisive role in the driving 



mechanism of pulsar glitchs, due to pinning of the neutron superfluid to the nuclei Jill]. Core 
nucleon superfluids may signiflcantly suppress cooling rates that rely on neutrino emission. 



by reducing the available phase space in the flnal state |12,13 
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In this work we focus on the inner core of neutron stars, where baryon species other 
than nucleons are expected to appear. It is widely accepted [14-18] that hyperons begin 
to accumulate at a density of about 2po, and at a density of 3po the hyperon fraction 
is already about 0.2. These results are a direct consequence of using modern estimates 
of the interactions of hyperons in nuclear matter, derived from hypernuclei experiments. 
The presence of hyperons has been shown to be of considerable importance in neutron star 
cooling rates due to their potential to participate in the efficient direct Urea processes. While 
hyperon direct Urea is found to be small compared to the nucleon direct Urea when nucleons 
are non-superfluid, hyperon direct Urea becomes the predominating coolant if the nucleons 
form superfluid pairs [0. It is noteworthy that the direct Urea mechanism can proceed 
through hyperon processes for almost any hyperon fraction, while the nucleon direct Urea 
requires a proton fraction of at least 0.11—0.15 [^,^]. In fact, some studies have found that 
hyperon direct Urea cooling is too rapid to be consistent with observed surface temperatures 
of pulsars [|l9i^^ . However, if hyperons also couple to a superfluid state, as expected for 
the nucleons, hyperon direct Urea will also be suppressed, and a large hyperon fraction could 
be easier to coincide with observed cooling rates. 

Hyperon pairing has not been studied previously, as the basic obstacles relevant to nu- 
cleonic pairing are pronounced for hyperons. However, a few measured events in KEK ex- 
periments attributed to doubly-strange AA hypernuclei, do offer indication with regard 



to the A — A interaction with a background nuclear matter medium. In a recent work Lan- 
skoy and Yamamoto formulated a G- matrix parameterization for the A — A interaction, 
based on Nijmegen OBE models. This G-matrix includes a dependence on the density of 
the nuclear matter medium, and reproduces the experimental results of the AA hypernuclei. 

In this work we aim to employ this formulation to estimate AA pairing energies in dense 
matter. We briefly review in Sec. II the formalism leading to the gap equation in the ^5*0 
channel. The properties of the effective potential used in this work are introduced in Sec. HI. 
Sec. IV presents our results for the superfluid gap of AA S-wave pairing in nuclear matter. 
Implications for neutron stars are discussed in Sec. V. Sec. VI contains our conclusions and 
some out-looks regarding hyperon pairing. 



II. THE GAP EQUATION 



The BCS theory predicts a transition to the superfluid phase when correlations 
leading to Cooper ^5*0 pairs give rise to excessive binding energy, which overcompensates the 
increase of energy due to the depopulation of the Fermi sea. The appropriate equations have 
been formulated in many works (see, for example, in refs. and for completeness we 
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review below the main results. We note in passing that variation of definitions may lead to 
differences in the numerical coefficients with respect to other works. 

The binding energy of a pair with momenta (k, — k) is found through a nonzero solution 
to the gap equation 



kk' 



2v {ii+^ir^ 

where is known as the gap function. The potential Vkk' is defined through the matrix 
element of the ^5*0 component of the interaction, and corresponds to the single particle 
energy, e^, when measured with respect to the Fermi surface. 

Going over to formal integration, the potential term is replaced by the potential matrix 
element: < k f — k [ \V\\i! } — k' |>. In the special case of the ^5*0 channel the matrix 
element is independent of the orientation of k and k'. For a two-particle central potential, 
V{r), the matrix element can be reduced to the form: 

Vuk' =< k\VCSo)\k' >= An / drUkr)Vir)joik'r) (2) 

Jo 

For convenience unit normalization volume is taken for the plane wave single particle wave 
functions; summation over spatial and spin exchange terms is implied. 

The integral form of the gap equation is thus: 

A, = -l-J-:^ I Ank'^ dk' Vkk'- ^" 



~ ^' ^^'vr^ TTT^ / ^ dr jo{kr)V{r)jo{k'r) . (3) 

In this work we use the common "decoupling approximation", where the Fermi surface 
is taken to be sharp even in the presence of of the pairing correlations. The functions are 
then simply given by 

Ck = £k — ^kp : (4) 



where we calculated the single particle energies with first order Hartree-Fock corrections 27 



The effect of the pairing potential on the single particle energies is often characterized 
by an effective particle mass, M*, which is typically lower than the initial (bare) mass by 
several percent. This mass can be estimated through the effective mass approximation: 
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which is usually found to be good up to a few percent 

Note that this effective mass differs from the bulk effective mass, found in field theories 



due to the meson scalar field, also typically lower than the bare mass [|T^]. A consistent theory 
thus requires an appropriate "true" initial mass, which includes medium effects through 
both the A— A and A— nucleon interactions. However, in the present work we invoke a non- 
relativistic approach, which has no means to consistently combine effective bulk masses, and 
correspondingly set the initial mass to be equal to the bare mass, i.e., M\ = 1115.6 MeV 
(some justification for this may also be found in uncertainties regarding values of effective 



masses at the Fermi surface 128|). The sensitivity of A pairing to this assumption is examined 
below. 



III. THE A - A POTENTIAL 



In this work we approximate the ^5*0 component of the A A interaction through a Brueck- 
ner G-matrix potential. We use the very recent G-matrix parameterization of Lanskoy and 
Yamamoto [^| derived from the Nijmegen OBE potentials for a AA pair in nuclear mat- 



ter. Their evaluation of the AA interaction is based on measurements of doubly strange 



hypernuclei observed in experiments p4| . Analysis of these experiments has suggested both 



the existence of an attractive component in the A — A interaction and the dependence of 



this interaction on the properties of the core nucleus The strength of the interaction 

is derived from the bond energy of the AA pair, defined as Ai?AA = -Baa — 25a- Here -Baa 
is the separation energy of two A's from the nucleus and -Ba is the separation energy of a 
single A from the same nucleus. 

The dependence of the interaction on the nuclear matter density is represented in Ref. 
[ p5|] by a three-range Gaussian form: 

3 

^AA(r) = + ^M^) + CikU^)) exp(-rV/3- ) , (6) 

i=l 

where kpin) is the nucleon Fermi momentum. Assuming symmetric nuclear matter (as is 
the case for light hypernuclei), kp is related to the nuclear density p^- by kF = (Stc'^^pn)^ ■ 
Since the A is an isospin singlet, it also seems safe to apply equation (^) to non-symmetric 
nuclear matter of density pn- 



The ranges (3i and the strength parameters ai,bi,Ci are taken from model ND of ||2^ 
and are listed in table |. This model successfully reproduces the experimental result of 
ABaa = 4.9 ± 0.7 MeV of 
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The radial dependence of the A — A interaction is demonstrated in Fig. 1 which shows 
VAA(r) for nuclear matter densities of Pn/po = 1, Pn/po = 2.5 and Pn/po = 5. At short 
distances the interaction is always repulsive, reflecting the core repulsion of the bare inter- 
action (we note that G-matrix approximations typically yield soft cores which substitute 
the need for short range cut-off necessary in other interaction models). At intermediate 
distances the ^So yields an attractive force of several tens of MeV's, which is strong enough 
to yield the pairing of the superfiuid state. 

As can be seen in Fig. 1, the dependence of the interaction on the nuclear matter density 
is rather weak. This implies that the existence of ^5*0 pairing should have only a moderate 
dependence on the density of the nuclear matter medium. We note, however, that the 
magnitude of the interaction tends to grow larger for a larger background density. 

It must be noted that the G-matrix parameterization is fitted to match experimental 
results for different nuclei, and is thus likely to be valid for a nuclear matter background 
with a density of pn^Po- In the following analysis we assume that the G-matrix is valid for 
higher densities as well. Clearly this is a somewhat crude assumption, especially since the 
G-matrix does not incorporate any relativistic effects which could be significant at densities 
relevant to neutron star cores (p>2po)- Hence, the results derived below must be viewed as 
preliminary estimates. More founded results must await better established AA potentials in 
high density nuclear matter. 

IV. RESULTS FOR A PAIRING IN DENSE MATTER 

Using the AA potential described in the previous Section, we have solved the gap equation, 
Eq. d^), for A hyperons in a nuclear matter background. The solution is found by iterations, 
when the integration is performed with a few hundred integration points, exponentially 
spaced around kpiA). The exponential spacing is required since the integrand in Eq. (H) 
is sharply peaked at the Fermi momentum. This behavior is demonstrated in Fig. 2 which 
shows the integrand for a A Fermi momenta of /cf(A) = 1.0 fm~^ and nuclear matter densities 
of p7v = 2.5po and pN = 5.0po. Fig. 2 also indicates the need for a large cutoff momenta in 
the calculation of the gap function in Eq. (^). 

Solution of the gap equations gives the gap function for any combination of values 
for the nuclear matter background density and the A Fermi momenta. Fig. 3 shows the gap 
function for the same values of pn and kpiA) as in Fig. 2. The gap function falls off from 
its maximum at k{A) = 0, and varies very rapidly around fcp'(A). The gap function is also 
found to be negative over a wide range of higher momenta. As can be seen in Fig. 3, the 
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gap energy is always larger in absolute magnitude for larger nuclear matter density. This 
results from the enhancement of the two-particle interaction at higher background densities, 
as seen in Fig. 1. Since the size of the superfluid gap for a given kplA) is mostly dependent 
on the two-particle interaction at distances of about l/kp{A), the gap energy grows larger 
for a larger density of the background nuclear medium. Correspondingly, The integrand of 
Eq. (0) (Fig. 2) also increases for a larger p^. 

We note that qualitative results such as those shown in Figures 2-3 are common also in 
solutions of the gap equations for nucleons The need for a large cutoff momenta is 

of particular importance, since it clarifies why the gap energy estimated through the weak- 
coupling-approximation (WCA) ||3l| systematically underestimates the gap energy. In this 
approximation one essentially assumes that it is sufficient to integrate Eq. (|^) over a narrow 
range near kp. Indeed, gap energies found for nucleons through the WCA are usually lower 
by a factor of two and more than those derived by a self consistent solution of the gap 
equation (see, for example, in Ref. |^). 

The prevalent result of the solution of the gap equation is the value for the gap energy 
at the Fermi surface, Ap = A^^. The resulting function Ap{kp) has a typical bell shape, 
ranging from Aq = to some maximum value and then falling off again to zero. This later 
part of the Ap{kp) arises from the decrease of mean inter-particle distance at higher kp, 
as the A's sample more of the repulsive core. This physical mechanism causes the S-wave 
superfluidity to vanish at large A partial densities. The Ap{kp{A)) dependence for nuclear 
matter background densities equal to 2po, 2.5po) 3po and 5po are shown in Fig. 4. The 
corresponding values of the gap energies and the effective A masses for = 2.5po and 
Pat = S.Opo are given in Table 2. 

As is expected from Figs. 2-3, the gap energy for a given kp{A) increases along with the 
density of the nuclear matter background, pn- However, for matter composed of nucleons and 
A's, increasing p^ alone corresponds to lowering the fraction of the A's. On the other hand 
keeping the A fraction constant while increasing the total density amounts to an increase of 
A Fermi momenta, and it is clear from Fig. 4 that increasing ^^(A) beyond 0.8 fm~^ should 
lead to a decline in the gap energy. Thus, increasing the total baryon density with a given 
A fraction tends to reduce the gap energy, while a larger total baryon density also means a 
larger nuclear matter density, which should increase the gap energy. Hence, these two trends 
compete when the total baryon fraction is increased and the A fraction is kept constant. 

In Fig. 5 we compare these two trends by presenting gap energies at the Fermi surface as a 
function of the total baryon density, ps, of matter composed of nucleons and A's. The curves 
represent constant A fractions of 5%, 10%, 15% and 20% of the total baryon population. As 
it happens for the particular pairing interaction used in this work, the two trends balance 
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for a A fraction of about 5%, and for a larger A fraction the gap energy decreases to zero 
as the total baryon density is increased. These results have direct implications on the gap 
energies in neutron star matter, as is discussed in the next Section. 

We now return to the problem of the "true" effective masses of the baryons in dense 
matter. So far we have assumed that the initial mass of the A hyperons on the Fermi surface 
is equal to the bare mass, Ma = 1115.6 MeV. Since we do not combine a self consistent 
treatment of bulk effects and the relativistic properties of the interactions, we must resort to 
arbitrary parameterization to examine the dependence of the pairing energies on the initial 
mass. A more accurate derivation of consistent interactions and masses is deferred to future 
work. 

Fig. 6 demonstrates the dependence of the gap energies on the "true" effective mass 
of the A hyperons in the matter. The results shown are for a nuclear matter density of 
PAr = 2.5po and A initial masses taken as 0.7, 0.85 and 1.0 times the bare mass (as mentioned 
in Sec. II, the effective mass derived by the solution to the superfluid equations is always 
lower than the initial one by several percent). As is expected, a lower mass leads to higher 
single-particle energies for any given momenta, and this yields lower gap energies. However, 
the basic existence of a superfluid gap of Af>0.1 MeV for /ci?(A) < 1.3 fm~^ is found also 
for effective masses lower than the bare mass. 

A final point of interest is the the dependence of the gap energies on the matter tempera- 
ture. The importance attributed to this dependence is of obvious in view of the implications 
of baryon superfluidity on neutron star cooling rates. We hereby follow the approach of 
Elgar0y et al. |^ in estimating this dependence for AA. 

The gap equation at a finite temperature T is given by revising Eq. (|^) to the form: 

A,(T)^-l' r 4-*'^ . tauh ( I , (7) 

2(2nfl (e(r) + A^,(T))i 1^ ikgT I' 

where ks is the Boltzmann constant. We solve Eq. (^, while approximating the single- 
particle energies to be "frozen", i.e., assuming that ^^/(T) =^fc/(0). This should be a reason- 
able approximation for neutron stars, since the temperature range of interest is much lower 
than the Fermi energy (see @] and references therein). We also assume that the two-particle 
interaction is not sensitive to the temperature in the range of interest. The gap equation is 
then solved in similar fashion to the zero-temperature case. 

The temperature dependence of the gap energy at the Fermi surface, Af{T) for back- 
ground nuclear matter density of piv = 2.5po and pN = 5po is shown in Fig. 7. Also shown 
are the critical temperatures, T^, estimated from the WCA, given by pl| 
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keTc ^ 0.57Af{T = 0) 



(8) 



As in the case of nuclear matter we see that the WCA does yield good agreement 
with the results of the full solution, provided that the value of Ap'(T = 0) is taken from the 
gap equation solution rather than the WCA for the gap, as explained above. 



V. IMPLICATIONS FOR NEUTRON STAR MATTER 



Modern estimates [14-18] of hyperon formation in neutron stars agree that hyperons 
begin to accumulate in neutron star matter at baryon densities of about 2po- In particular, 
the threshold baryon density for A formation is found to be about 2.5po, when the chemical 
potential of the neutrons grows large enough to compensate for the mass difference Ma — M„. 
While the fine details of the A fraction in the matter are model dependent, these basic features 
are widely accepted. We stress that this consensus is an immediate result of employing 
realistic values for the interaction of A hyperons in nuclear matter, based on experimental 
data of A-hypernuclei 



An example of the equilibrium composition of neutron star matter (assuming T = 0) is 
given in Fig. 8a, based on an equation of state similar to the ^ = 7 = | model of Ref. |jl8|. 
The steep rise in the A fraction when they first appear in the matter is common to all works 
that examined hyperon formation in neutron stars. This behavior is caused by the fact that 
lowering the nucleon fraction lowers the nucleon-nucleon repulsion and the nucleon Fermi 
energies, while the net interaction among the A's is still attractive. Eventually the A fraction 
saturates, typically at 0.1 — 0.2, and continues to grow slowly up to as much as 0.3 at higher 
densities. 

Recent theoretical and experimental results of S^-atoms suggest that the interaction of 



S hyperons in nuclear matter includes a strong isoscalar repulsive component ||32|. If such 
repulsion exists, formation of S hyperons in neutron star matter is suppressed [|16|,|18[, and 
A production in the matter is somewhat enhanced, both by a lower threshold density and by 
a sharper rise of the A fraction. The main effect, though, is the formation of S~ hyperons 
which begins at significantly lower densities (about 3po), providing the favorable negatively 
charged baryon fraction. The equilibrium compositions of matter without S's is shown in 
Fig. 8b, using an equation of state otherwise identical to that of Fig. 8a. 

In view of the absence of any experimental data on medium effects regarding different 
hyperon species, we assume in the following analysis that the pairing interaction discussed 
in Sec. 3 is valid also for a background matter which includes other species besides nucleons 
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(i.e. S and S hyperons). For densities up to ~5po this is a reasonable assumption, since the 
non-A matter is highly dominated by the nucleons. Thus, for every combination of the total 
baryon density and particle fractions we take Pbg = Pb—Pa as the background density pjv for 
the calculation of the ^5*0 gap energy. 

The AA gap energies found for the baryon compositions of Figs. 8a-8b are shown in 
Fig. 9, as a function of the total baryon density. Also shown are the gap energies for the 



equilibrium composition for model PLZ of Schaffner and Mishustin |17], which predict A 



accumulation at slightly lower densities than the equations of |18 



As seen in Fig. 9, the qualitative behavior of the gap energies is common to all three 
equations of state. Pairing to a superfluid state essentially takes place once the A's appear in 
the matter, and rises sharply to a maximum value following the sharp rise of the A fraction 
in the matter. The partial density of the non-A baryons is almost constant in this range of 
total baryon densities. Hence, the curves approximately follows the gap energy dependence 
on kpiA) for a given background density, as shown in Fig. 4. The pairing energy rises 
sharply as the total baryon density is increased, and then, as is expected from Fig. 5, begins 
to decline once the A fraction exceeds about 0.05, (note that for pB = 2pQ, kpiA) ^0.8 fm~^ 
is reached when the A fraction is about 5%). Since the A fraction begins to saturate at a 
value of 0.1—0.2, the decline of the gap energy is not as steep as in the rising part. The 
rate of this decline is thus somewhat model dependent, particularly whether other hyperon 
species (i.e. the S~) compete with A formation. 

It should be noted that these results are qualitatively similar to those found for proton 
^5*0 pairing in neutron star matter, where protons are a minority among the nucleons. The 
density range found for a superconducting proton state lies between the threshold for free 
proton appearance up to densities where the proton fraction reaches about 0.1 — 0.2 
However, since the proton fraction in neutron star matter is expected to rise much more 
moderately as a function of the total baryon density than the A fraction (see Figs. 8a-8b), 
the density range where a proton superconductor exists is typically larger than that found 
here for the A superfluid. 



VI. CONCLUSIONS AND OUTLOOK 



In this work the ^5*0 pairing energy of A hyperons in a nuclear matter background was 
evaluated using the G-Matrix effective interaction presented by Lanskoy and Yamamoto ||25|| . 
We find that a gap energy of a few tenths of MeV is expected for a A Fermi momenta, /cir(A), 
below 1.3 fm~^. The gap energy is dependent both on the A Fermi momenta and on the 
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density of the background nuclear matter, p^. For pat > 2po the gap energy for a given kp^^s) 
increases with increasing p^. 

Employing these results to neutron star matter with hyperons yields A A ^5*0 pairing 
for a baryon density range between the threshold density for A appearance to about the 
baryon density where the A fraction reaches ~ 0.2. A maximum gap energy of 0.8 — 0.9 
MeV is achieved for a A fraction of about 0.05. While the exact range of densities where 
such pairing exists is model-dependent, the qualitative picture seems to be common to all 
equations of state which are based on modern evaluations of the A— nucleon interaction in 
nuclear matter. Gap energies in this range are larger than the temperature predicted in 
neutron star cores, and thus imply that a A ^5*0 superfiuid will exist in the core, typically 
within a baryon density range of pb~2 — 3po- 

We comment that the present results must be treated as a preliminary evaluation of A 
pairing in dense matter. The evaluation of the two-particle interaction is based on hypernu- 
clei experiments, where the nuclear matter density is limited to Pn^Po, so that the effective 
interaction might not be as good an approximation as in the case of neutron pairing. In 
particular, the present work does not include relativistic corrections which might be signifi- 
cant at the baryon densities where hyperons form in neutron star matter (note however, that 
relativistic corrections for proton ^5*0 pairing at about the same densities have been found 
to introduce only small corrections to the nonrelativistic results [^). It is also noteworthy 
that we have not included particle-hole correlations which have been shown to be important 
in the evaluation of gap energies 0. In short, further work is necessary to produce more 
realistic results, preferably with a better founded A— A interaction in a high density nuclear 
matter background. 

We believe that formal treatment of the non-nucleon component in the background neu- 
tron star matter, will not significantly effect the results found here. This is especially true 
if S hyperon formation is suppressed, so that the baryon equilibrium compositions include 
only nucleons and A's throughout the entire range where pairing is expected. Nonetheless, 
taking other hyperon species into account is clearly desirable in a more rigorous model. Ob- 
viously, hyperon-hyperon interactions in a dense matter background will also provide a basis 
for estimation of possible pairing of other hyperon species in neutron stars. For example, 

pairing is of special interest, since the S~ is also expected to appear at relatively low 
baryon densities in neutron stars (if S formation is not suppressed). However, no relevant 
experimental data is currently available. The commonly assumed universal hyperon-hyperon 
interaction implies that the AA gaps may serve as indication for S and S pairing in dense 
matter. More accurate results require, however, the inclusion of isospin-dependent forces, 
which are absent in the A case. 
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The large majority of dense matter equations of state require neutron star central den- 
sities larger than the threshold density for A formation, i.e. baryon densities larger than 
~2.5po- Hence, it is likely that neutron stars do include a region where the A's pair to a ^5*0 
superfluid. Whether or not the central density of a neutron star exceeds the density range for 
^5*0 AA pairing depends on its mass and on the actual equation of state. Note, however, that 
at larger densities higher order pairing may also be available, including inter-species pairing 
[ p^ . In fact. An pairing may be more likely than pn pairing, since at baryon densities of 
Pb>^4:Po the A and neutron fractions are expected to be comparable. 

Finally, we recall that the existence of a ^5*0 A superfluid for baryon densities relevant to 
neutron stars implies significant suppression of A-direct Urea cooling. The onset of superflu- 
idity reduces the neutrino emissivity, along with the heat capacity and thermal conductivity, 
by a factor of exp(— A^/Aj^T). In view of our results here, we suggest that implications of 
hyperon superfluidity on neutron star cooling rates are well worth examination. 
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Figure Captions 

Figure 1: The radial dependence derived for the A A G- matrix interaction presented in 
[^. The curves correspond to nuclear matter background densities of pa? = 2.5po and 
5po, where po is the nuclear saturation density. 

Figure 2: The integrand of the gap equation, Eq. (0), for A;i?(A) = 1 fm^^, as a function of 
the secondary momenta k'. The curves correspond to nuclear matter background densities 
of pAr = 2.5po and Spo- 

Figure 3: The gap function, A^, when kpiA.) = 1 fm^^ and nuclear matter background 
densities of piv = 2.5po and Spo- 

Figure 4: The gap energy Ap for AA pairing as a function of the Fermi momenta, for 
nuclear matter background densities of pAr = 2po, 2.5po, 3po and 5po. 

Figure 5: The gap energy for A A pairing as a function of the total baryon density, pb, 
for different fixed A fractions. 

Figure 6: The gap energy Ap for AA pairing as a function of the Fermi momenta, 
for different values of the initial mass of A hyperons in the matter. The nuclear matter 
background density is taken as pAr = 2.5po. 

Figure 7: Temperature dependence of the gap energy for ^^^(A) = 1 fm~^ for nuclear 
matter background densities of = 2.5po and 5po. Also indicated are the corresponding 
weak-coupling estimates for the critical temperatures. 

Figure 8a: The equilibrium compositions of neutron star matter with hyperons, as a 
function of the total baryon density, p^. The compositions were calculated with an equation 
of state similar to the 5 = 7= | model from flsH . 

Figure 8b: Same as Fig. 8a, but when S hyperons are repelled by the nucleons and their 
formation is thus suppressed. 

Figure 9: The gap energy of A A ^5*0 pairing in neutron star matter as a function of 
the total baryon density. The equilibrium compositions of the matter are those of Figs. 8a 
(BG+E) and 8b (BG-S), and for model PLZ of Schaffner and Mishustin (SMpLz)- 
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TABLE I. 

Parameters of the ^5*0 state of the AA G-matrix potential (model ND of [p5|] ) 



A 






bi 




Ci 


(fm) 


(MeV) 




(MeV fm) 




(MeV fm2) 


0.5 


835.5 




-252.7 




122.7 


0.9 


-298.5 




156.6 




—55.07 


1.5 


-10.80 




3.0398 




-1.126 




TABLE IL 








A A ^5*0 pairing 


energy gaps and effective 


masses 










PN = 2.5po 








PN = 5P0 




M*/M 


Af 




M*/M 


Af 


(fm-i) 




(MeV) 






(MeV) 


0.2 


0.9967 


0.0432 




0.9963 


0.1321 


0.3 


0.9895 


0.1767 




0.9881 


0.3897 


0.4 


0.9771 


0.3749 




0.9740 


0.7143 


0.5 


0.9596 


0.5868 




0.9543 


1.0371 


0.6 


0.9383 


0.7628 




0.9304 


1.2998 


0.7 


0.9150 


0.8677 




0.9042 


1.4611 


0.8 


0.8915 


0.8735 




0.8779 


1.4937 


0.9 


0.8693 


0.7826 




0.8532 


1.3876 


1.0 


0.8497 


0.6130 




0.8313 


1.1574 


1.1 


0.8335 


0.4027 




0.8132 


0.8409 


1.2 


0.8211 


0.2053 




0.7993 


0.5262 


1.3 


0.8128 


0.0495 




0.7900 


0.1810 
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figure 6 
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